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Abstract 

We examine the properties of cold ions confined by a Paul trap in a linear crystal con- 
figuration, a system of considerable current interest due to its application to practical 
quantum computation. Using a combination of theoretical and numerical calculation, 
a semi-empirical formula for the positions of the ions is derived. 
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Cold ions confined in electromagnetic traps are rapidly becoming a very impor- 
tant system both for the study of fundamental physical systems, such as cold charged 
plasmas or quantum chaos, and for technological applications such as optical frequency 
standards. Recently a chain of cold ions in a linear trap has been proposed as a pos- 
sible means to realize a quantum computer |]J. This idea was confirmed in principal 
soon after when an elementary quantum logic gate was realized experimentally using 
a trapped Beryllium ion 0, and larger scale devices are currently being pursued by 
several experimental groups (see, for example, ||). Understanding the properties of 
collections of confined ions is of great importance to these endeavors. 

As is well known [f|], it is impossible to confine charged particles by electrostatic 
forces alone. To overcome this problem, the radio-frequency Paul trap was developed: 
such devices use a electromagnetic field varying at radio frequencies (~ 100 MHz say) 
to produce an effective binding potential in three dimensions [H. When two or more 
ions are confined in such a trap, they will repel each other due to the Coulomb force. As 
a result, such confined charged plasmas will have very low densities. When sufficiently 
cold, the plasma will condense into a crystalline state. In the highly anisotropic traps 
used for some atomic clocks and for quantum computing, this crystalline state is, 
for small enough numbers of ions, a simple chain of ions lying in a straight line. As 
the degree of anisotropy is decreased, or number of ions is increased, phase changes 
to other configurations will occur: firstly the ions adopt a zig-zag configuration, and 
then a helical configuration. These phase changes have been studies numerically 
and analytically ||. 

In this letter we present results of new numerical studies of the positions of ions 
confined in highly anisotropic traps. This information is of course of considerable 
importance in designing trapped ion quantum computers. Using a simple theoretical 
argument we then develop a relatively compact expression for the position of each ion, 
which depends on the total number of ions confined in the chain. Our results are 
compared with both our numerical data, and with results obtained previously by other 
authors, and good agreement is obtained. 

Consider a chain of N ions confined in a linear trap (fig.l). The position of 
the n th ion, where the ions are numbered from left to right, will be denoted by the 
position vector relative to the trap center (i.e. the minimum of the binding potential) 
R-n(^) = (A„(t), Y n (t), Z n {t)). The motion of each ion will be influenced by an overall 
harmonic potential due to the trap electrodes, and by the Coulomb force exerted by 
all of the other ions. Thus the potential energy of the ions in the ion chain is given by 
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the following expression 



M N 

F(R 1 ,R 2 ,...R 7V ) = f-E(^n + -X 2 + ^ 

z n=l 

8^6(1 , Rm — Rm 

u n,m=l 1 " " l| 

where M is the mass of each ion, e is the electron charge (the ions are assumed to 
be singly ionized), eo is the permittivity of free space and uj x is the angular frequency 
characterizing the strength of the trapping potential in the x direction (and similarly 
ojy and w x for the y and z directions). 

The equilibrium positions of the ions, R^ are defined by solutions of the following 
equations 

[W(Ri,R2,...Rar)] R _ R (0 = (2) 
Substituting from (|l|) we obtain 

^-^sW^§p =o(i=i,2,3) " (3) 

where i = (1, 2, 3) denote the X, Y and Z components, respectively. 

Let us assume that the trap potentials are sufficiently strong in the Y and Z 
directions and sufficiently weak in the X direction that in equilibrium the ions lie in a 
straight line along the X-axis. Mathematically this assumption is expressed by 

Ri°) = £( Mn ,0,0) (4) 

1 /3 

where £ is a scale length given by (e 2 /47re Afu;2) 1/J an d u n is the dimensionless equi- 
librium position of the n-th ion, which is a solution of the following set of N coupled 
algebraic equations, obtained by substitution from eq(f|) into eq(|3]): 

u m +J2 T iU - " V ] = (™ = 1,2,...N), (5) 

n= l \U"m U n ) 

where sgn(x) = 1 if x > and — 1 if x < 0. For N = 2 and N = 3 these equations may 
be solved analytically 0: 

N = 2: Ul = -(1/2) 2 / 3 , M2 = (l/2) 2 / 3 , (6) 
N = 3 : mi = -(5/4) 1 / 3 , u 2 = 0, u 3 = (5/4) 1 / 3 . (7) 
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For iV > 2 it is necessary to solve for the values of u m numerically. For small numbers 
of ions (N < 100 say) a Newton- Raphson method can be employed to find u m ; however 
this becomes inefficient as N gets large. Therefore we used another method, based on 
the following set of equations of motion 



><„, IT) = - , r . - v m[ r) + £ s g n M r ) ^(r)) „ ( m L; 2 V). (8) 

n=l 0w) ~ V n (T)) 

where the single and double dots denote single and double differentiation with respect 
to the dimensionless time variable r = u x t. These equations represent a hypothetical 
damped oscillation of the ions in the trap, including their mutual Coulomb interaction. 
The solutions of these equations have the property that 

lim v m (r) = u m , (9) 

T— >00 

where u m are the desired solutions of eq(||). The integration of eq® was carried 
out numerically using the standard fourth-order Runge-Kutta method|l(J. Because 
eq(^|) is valid regardless of the initial conditions, the simplest possible initial conditions 
were used, i.e. all of the ions being equally spaced. The Runge-Kutta algorithm was 
applied repeatedly until the values of u m between adjacent iterations were identical to 
the seventh decimal place. This was done for up to 1000 ions (although not for all 
numbers). The values there by obtained are in agreement with those obtained (for 
N < 100) by Newton- Raphson ||. This dynamic technique can be adapted quite easily 
to study classical wave motion in the ion chain; this will be the subject of a forthcoming 
paper. 

In order to make some sense of the large amount of data generated we will now 
derive a analytic formula which approximates the numerical results quite closely. Our 



analysis is based on the very elegant idea due to Garg [p| . Let us consider the force 
acting on an ion at position X. The Coulomb force due to the two nearest neighbor 
ions is 



F 

- 1 nn 




(10) 



4vre 5(X) 2 dX ' 

2 A data file called ion_positions.dat which contains the results of these numerical calculations 
can be found in the directory pub/james/Ion_Position_Data which can be accessed via anonymous 
ftp to t4 . lanl . gov. 
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where S- is the distance from the ion to the nearest neighbor on the left, S + is the 
distance to the nearest neighbor on the right and S(X) is the separation of ions at 
position X, treated as a continuous function. This is a reasonable approximation 
to make for large numbers of ions. The next nearest neighbors are approximately 
twice as far away as the nearest neighbors, and so the force they exert on the ion is 
approximately _F nn /4; the next pair of ions are three times as far away as the nearest 
neighbors, and so the force they exert is approximately F nn /9, and so on. Thus the 
total Coulomb force on the ion will approximately be given by the following expression: 

OO 1 

F C ~ fnnEp 
k=l K 

e 2 tt 2 dS(X) 



Ane 3S(X) 2 dX ' 

where we have used the fact that Y^kLi l/^ 2 — vr 2 /6 and we have approximated the 
finite sum over all ions as an infinite sum. This approximation should be valid near the 
center of the ion chain, but will not yield very good results at the ends of the chain. 

The Coulomb force acting on the ion at position X will be balanced by the har- 
monic restoring force due to the trap electrodes. Thus we can write the following 
identity: 

P 2 TT 2 J.Q(Y\ 

Mu 2 x X = 0. (12) 



47re 35(X) 2 dX 

If we introduce the dimensionless ion separation a(u) = S(X)/£ and dimensionless 

1 /3 

distance from the trap center u = X/£, where, as before t = (e 2 /47reoMa; 2 ) , we 
obtain the following differential equation for the separation: 

— = —ua . (13) 

dU 7T Z 

This can be solved quite easily, yielding the formula [ |TT|] : 



, , 27r 73 

= ^ (14) 



where C is a constant, which could be determined from the value of the separation of 
ions at the trap center (u — 0). 

Let n(u) be the total number of ions which are within a scaled distance u of the 
trap center. Clearly n(u) is given by the following differential equation 

du ^{u) 
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On substitution from eq(PD, and performing the integration, we obtain the following 
formula for n(u): 

n(u) = Au- Bu 3 , (16) 

where we have set n(0) = 0. The constants A and B can be related to the constant 
C introduced above. However, instead of attempting to carry this analysis too far, 
it is better at this stage to obtain empirical formulas for the the constants A and B 
based on our numerical results. This was done by performing a least squares fit of 
the numerical data to a cubic formula of the type given by eqflTB"!). The values of A 
and B were found for a variety of different total numbers of ions. When this data was 
compiled, we found that A and B were approximately given by the following power 
laws: 



A(N) w 0.436iV a596 
B(N) « 0.0375iV" a178 , 



where N is the total number of ions in the chain. 

To obtain an expression for the position of the n-th ion in the trap, it is necessary 
to invert eq fllB"!) . This can be done using the standard formulas for the roots of a cubic 
equation |I2[ . We therefore obtain, taking care to select the correct root based on the 
value of n at u = 0, the following formula for the the scaled equilibrium positions of 
the n-th ion: 



u, 



UA 
W 



cos 



COS 



I27B 
4^3 



n — 



(N+1) 



+ 



Air 



a(N) sin ( - sin" 



P(N){n- 



(N+l) 



(1* 



If we reintroduce the scale length £, we finally obtain the following expression for the 
equilibrium position of the n — th ion, when there are a total of N ions in the trap: 



x (o) 



1/3 



a(N) sin f - sin" 



P(N)\n 




(19) 



where, as before, the ions are numbered from left to right, M is the mass of each 
ion, e is the electron charge, €q is the permittivity of free space and uj x is the angular 
frequency characterizing the strength of the trapping potential in the x direction. In 
eqsflg) and (0) we have introduced the coefficients a{N) = ^4A/3B w 3.94iV - 387 
and P(N) = J 27 B/AA 3 w 1.75A^-°- 982 . 
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Equation ([19]) is the main result of this note. As an example we have plotted in 
Figure 2 the numerically calculated ion positions together with the positions calculated 
using this formula, for a total of 41 ions in the trap. Also we have included experimental 
ion position data gleaned from Fig.5 of reference |§. As can be seen from the figure, 
there is good agreement between the numerical data and the empirical formula. The 
differences between the experimental data and that calculated numerically, may well 
be due to the departure of the trapping potential from the harmonic form we have 
assumed. The percentage r.m.s. error between the ion positions calculated numerically 
and those calculated using eqfll~8|) is shown in Figure 3. The error is only of the order 
of a few percent when N > 25, but as expected, errors increase for small numbers of 
ions. 

For small arguments one can make the approximation sin(sin _1 (x)/3) ~ x/3, and 
so, near the trap center [where (n — N + 1/2) is a small number] the scaled ion positions 
are given by: 

u n « 2.29 (n-N+ 1/2) jy- - 596 . (20) 

Hence the minimum separation between ions, which is of considerable importance in 
quantum computer design is given by 

u min (N) « 2.29iV- a596 . (21) 

This result is in good agreement with the empirical formulas previously calculated for 
the minimum separation of ions are the trap center [13], [14]|. This formula is plotted 



in figure 4, along with the numerical data, and the following formula for the minimum 



ion spacing due to Dub in gj (see also 0]), based on a fluid model for the ion cloud: 



t (N) « 1.92iV- 2/3 ln(aiV) 1/3 , (22) 



where a = 6e 7 " 13/5 « 0.794, 7 being Euler's constant. As can been seen from figure 
4, both the empirical formula derived here and the analytic formula due to Dubin 
approximate the numerical data quite closely. 
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Figure Captions 

Figure 1. A schematic illustration of ions confined in an harmonic trapping potential. 

Figure 2. Equilibrium positions of ions when there is a total of 41 in the chain, as 
calculated numerically (crosses) and by eq (|i~8|) (plane line). Also show are the exper- 
imental positions of Hg+ ions in a linear trap, which were gleaned from fig.5 of || 
(circles). 

Figure 3. Root mean square percentage error for calculating positions of trapped ions 
using formula eq(|T8|) for total ion numbers N up to 1000. 

Figure 4. Comparison of numerical results for minimum ion separations (crosses) with 
the empirical power law eq([H|) (plane line) and the analytic formula eq (|22|) (dashed 
line). The two curves are in such good agreement that it is difficult to distinguish 
them. 
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